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The solar wind wake behind the moon is studied with ID electrostatic particle-in-cell (PIC) simulations using 
a physical ion to electron mass ratio (unlike prior investigations); the simulations also apply more generally 
to supersonic flow of dense magnetized plasma past non-magnetic objects. A hybrid electrostatic Boltzmann 
electron treatment is first used to investigate the ion stability in the absence of kinetic electron effects, 
showing that the ions are two-stream unstable for downstream wake distances (in lunar radii) greater than 
about three times the solar wind Mach number. Simulations with PIC electrons are then used to show that 
kinetic electron effects can lead to disruption of the ion beams at least three times closer to the moon than 
in the hybrid simulations. This disruption occurs as the result of a novel wake phenomenon: the non-linear 
growth of electron holes spawned from a narrow dimple in the electron velocity distribution. Most of the 
holes arising from the dimple are small and quickly leave the wake, approximately following the unperturbed 
electron phase-space trajectories, but some holes originating near the center of the wake remain and grow 
large enough to trigger disruption of the ion beams. Non-linear kinetic-electron effects are therefore essential 
to a comprehensive understanding of the ID electrostatic stability of such wakes, and possible observational 
signatures in ARTEMIS data from the lunar wake are discussed. 

PACS numbers: 52.30.-q, 52.35.Fp, 52.35.Mw, 52.35.Qz, 52.35.Qz, 52.35.Sb, 52.35.Tc, 52.65.Rr, 95.30.Qd, 
96.20.-n, 96.25.Qr, 96.25.St, 96.50.Ci 


I. INTRODUCTION 

The supersonic flow of the solar wind past the moon 
and other objects without global magnetospheres is a 
topic of continued theoretical and experimental interest. 
The ongoing ARTEMIS missiorP is studying the solar 
wind in the vicinity of the moon, with targeted 3D hy¬ 
brid simulations being usecP to interpret the observed 
fluctuation^^. More generally, simulations of dense mag¬ 
netized plasma flow past non-magnetic objects has wide 
applicability, ranging from probes in laboratory plasmas 
to asteroids and the moons of SaturrP^. 

Solar wind flow past the moon leads to an elaborate lu¬ 
nar plasma environment and wak^. At distances beyond 
a few lunar radii the main dynamics of interest are those 
associated with the ambipolar plasma expansion into the 
void trailing the moon, and in particular any resulting 
instabilities. This process has been studied with fully 
kinetic simulations in IECHU 

and 2rPi, a s we ll as with 
hybrid simulations in IDEI, 2EP^, and 3 d[ 1II3HII1, There 
are also 2D kinetic si mulati ons which ignore the magnetic 
field of the solar wind^^^, but these are less relevant to 
the moon given its large size relative to the gyroradius 
scales in the solar wind. 

Despite the wide range of prior simulations, none have 
so far accurately captured the ID electrostatic phenom¬ 
ena which can occur in the wakes of large non-magnetic 
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objects. Birch and ChapmaiP attributed the fluctuations 
they (and previously Farrell et al.^^ observed to elec¬ 
tron beam instability. However, HutchinsorP^ showed 
that the unphysical ion to electron mass ratio used in 
those simulations artificially enhances the energy scale 
of electron-electron instability, indicating that a realistic 
mass ratio is needed to accurately model such instabilities 
in the wake. The present work confirms the importance 
of using a realistic mass ratio, discovering a novel wake 
phenomenon where electron holes grow from the small en¬ 
ergy scale of the electron-electron instabilities to a scale 
where they can disrupt the counterstreaming ion beams 
in the wake. This observation has already inspired the¬ 
oretical investigations into the underlying mechanisnP^, 
which show that the changing density in the wake can 
act as a drive for hole growth. 

Kinetic electron effects are thus shown to be impor¬ 
tant in the flow of magnetized plasmas past objects much 
larger than the plasma Debye length. The case of super¬ 
sonic flow is treated here, but the failure of Boltzmann 
or other hybrid electron treatments to capture the rele¬ 
vant electrostatic phenomena may well extend to cases 
with slower plasma flow: counterstreaming ions will also 
be present for slower flow, and the electrons will still see 
a varying potential and local density, so the main dif¬ 
ferences are geometrical and related to boundary condi¬ 
tions at the collecting object. Whether these differences 
are enough to render the kinetic electron effects unim¬ 
portant is a topic that should be investigated, as it may 
be relevant to Mach probes and other means of invasive 
plasma measurements. 

This paper is organized as follows: The ID electro¬ 
static simulation approach used is described in Section[n| 
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Simulations using Boltzmann electrons are presented in 
Section m showing the onset of ion two-stream insta¬ 
bility; the unstable (but mass-ratio dependent) electron 
distributions associated with the Boltzmann-simulated 
wake are also illustrated. Fully kinetic simulations are 
presented in Section [IV| where some of the electron holes 
formed early in the wake are seen to grow and lead to dis¬ 
ruption of the ion beams much closer to the moon than 
in the Boltzmann-electron simulations. Possible observa¬ 
tional signatures of electron holes and ion-beam disrup¬ 
tion in ARTEMIS data are discussed in Section |V| and 


some concluding remarks given in Section VI 


II. ID ELECTROSTATIC SIMULATION METHOD 

The lunar wake simulation is reduced to ID through 
the same approach taken in prior such studieJ^® by fol¬ 
lowing the plasma in the frame moving with the solar 
wind, where the perpendicular motion is restricted by 
the magnetic field. What is simulated is then the par¬ 
allel motion of electrons and ions along an unperturbed 
magnetic field line as a function of time. For a static 
wake this time-dependence can be translated to a spa¬ 
tial dependence by following the trajectory of the field 
line, and a single ID simulation thus applies to a wide 
range of solar wind flow velocities and magnetic field ori¬ 
entations with respect to the moorP^ (in suitably scaled 
coordinates). 

For ease of comparison with the prior studie^^K^ the 
relative motion of the solar wind and the moon is taken 
to be perpendicular to the magnetic field and have speed 
'T’sw = 25cs. Here Cg = is the (cold ion) sound 

speed, which is equal to the ion thermal velocity since 
Te = Ti. The distance behind the moon is thus x = 
and the position parallel to the magnetic field is denoted 
y. The particular choice of solar wind speed does not en¬ 
ter into the simulations, only the representation of the re¬ 
sults in the figures, so other speeds can be accommodated 
by a simple scaling of x. Note that the convention for x 
used herein differs from that used in References M and 
HHby a factor of the (cold ion) Mach number ^sw/cg = 25. 

A ID electrostatic particle-in-cell (PIC) code, referred 
to here as ESPIC, has been developed to carry out the 
simulations. The ions and electrons are evolved with a 
standard leap-frog scheme at time steps small enough to 
resolve the electron motion, and their charge distributed 
to the two adjacent nodes of the mesh using the cloud- 
in-cell approaclP^. At each time-step the finite-difference 
Poisson equation is solved for the potential using a direct 
tridiagonal method, yielding a constant electric field be¬ 
tween any two mesh nodes to accelerate the particles. 
The potential at the boundaries is set to zero with ho¬ 
mogeneous Dirichlet boundary conditions, and particles 
are injected there according to the assumed Maxwellian 
background distributions. 

ESPIC can be run in parallel to accommodate the large 
number of particles and time-steps needed to resolve the 


wide ranges of relevant length and time scales. Even so, 
large scale-separations such as in the case of the lunar 
wake present a challenge: the lunar radius Rm is 1730 km, 
while the ion thermal gyroradius is around 40 k m and the 
electron and ion Debye lengths are around 7 npIII2|_ 
resolve the Debye length in a domain many lunar radii 
across would thus require ^10^ grid-points, allowing for 
a few grid-points per Debye length to reduce discretiza¬ 
tion error. That is not currently feasible using a realistic 
ion to electron mass ratio, but extrapolation based on 
simulations at larger (but still small) Debye lengths can 
yield insight. 

For such ID electrostatic simulations to be applicable, 
the main influence of the magnetic field on the particles 
must be to restrict their motion in the directions per¬ 
pendicular to it. Further, the collecting object (e.g., the 
moon) must have small apparent extent in the direction 
of flow, such that the ions do not move much as the field 
line traverses the object. This can either be because the 
flow is sufficiently supersonic or because the object has a 
large aspect ratio, and in either case the detailed shape of 
the object doesn’t matter; its main influence is to leave a 
region depleted of plasma behind it. The parallel dynam¬ 
ics of the ambipolar plasma expansion into that depleted 
region then dominate the solution, and can be captured 
in a ID domain by suddenly removing the ions and elec¬ 
trons from the part of that domain crossed by the object; 
this is the approach taken in ESPIC. 

Despite not directly influencing the ions, the shape 
(and surface charge) of the object could still affect the 
(much faster) electrons. For example, in the case of the 
lunar wake the electron thermal speed is faster than the 
solar wind flow (for a realistic mass ratio), so it is not 
very accurate to model the passage of the moon by a 
sudden removal of electrons. However, the simultane¬ 
ous removal of both species provides a convenient way 
to remove an equal number of electrons and ions, which 
avoids unphysical charge build-up associated with inac¬ 
curate object boundary conditions (especially important 
for simulations with short Debye length). Further, as 
will be shown in Section m the initial electron dis¬ 
turbance (and thus any inaccuracy associated with it) 
quickly propagates out of the domain in simulations with 
a realistic mass ratio. Thus, the way of implementing the 
(2D) passage of the moon in the ID ESPIC simulations 
is thought to accurately capture the relevant effects on 
the part of the wake of interest in the present work. 

The main differences between the ESPIC simulations 
and prior ID kinetic simulationJ^t^ are finer phase-space 
resolution (up to 2 x 10^ grid points and 10^ particles), 
the use of a physical mass ratio, and the absence of pe¬ 
riodic boundary conditions. Preventing the wake from 
interacting with itself through the periodic boundaries 
set a lower limit on the domain size in the prior simula¬ 
tions, and for a physical mass ratio that constraint would 
be unmanageable. Care must still be taken to avoid re¬ 
flection of potential structures from the boundaries, but 
that is much less restrictive in terms of domain size than 
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the periodic boundaries; if need be absorbing potential 
boundary conditions could be added allow even smaller 
domains. 

For small Debye length Ad, a random initialization of 
particles in the PIC method yields unphysically large 
fluctuations in the initial potential. This is especially 
true at long wave-lengths, as the amplitude of potential 
fluctuations of wavelength A can be shown to scale with 
for Gaussian counting statistics. For the same 
number of particles, the amplitude of the potential fluc¬ 
tuations can be dramatically reduced by using a quiet 
starW^^ where particle positions and velocities are ini¬ 
tialized more uniformly than in a true random start; a 
quiet start is used in ESPIC. 


III. BOLTZMANN ELECTRON SIMULATIONS 


Hybrid simulation approaches involve the use fluid or 
other non-kinetic treatments of electrons to greatly re¬ 
duce the computational cost of simulations, which then 
only need to resolve the ion time-scales as opposed to 
the (typically much shorter) electron ones. One such ap¬ 
proach is to assume that the electron density along a 
magnetic field line satisfies a Boltzmann relation 


rie = neo exp 



( 1 ) 


which (for instance) it would exactly along a static, 
monotonically decreasing (repulsive) potential (j) from 
some point with a stationary Maxwellian electron dis¬ 
tribution of density neo- Such Boltzmann electrons are 
implemented in ESPIC by solving the non-linear Poisson 
equation arising from using Equation for the electron 
density, and are used in the simulations throughout this 
section (even when moving electrons kinetically to exam¬ 
ine their distribution). 

The wake structure resulting from a simulation with 
Boltzmann electrons and a solar wind electron Debye 
length Aoe = 0.02Rm is shown in Eigure[2 As discussed 
in Section [ll] a solar wind speed Ugw = 25cs is used to fa¬ 
cilitate comparison with prior ID kinetic simulation^^t^. 
The Boltzmann-electron simulation captures the main 
features of the initial plasma expansion into the void 
behind the moon seen in the prior simulations: A neg¬ 
ative potential in the density-depleted region approxi¬ 
mately balances the electron and ion densities, with sig¬ 
nificant charge density seen only in the transition from 
unperturbed plasma to the strongly density-depleted re¬ 
gion. Later in the wake, for x ^ 140i?M, the simulation 
shows rapid ion (and charge) density variations (verti¬ 
cally elongated striations in the plots) resulting from ion 
two-stream instability ( somet imes referred to as the ion- 
ion [acoustic] instabilit}®^ . 

Eigurej^ shows the ion phase space for the Boltzmann- 
electron simulation at three different times/positions. 
The main initial response of the ions is that they are 


accelerated towards the void by the negative potential 
structure. Later in the wake free-streaming of ions is 
the dominant process, distorting the initial phase space 
void such that the density-depletion is shallower but more 
spread out. This sets up a configuration of two ion beams 
across most of the wake, and the velocity gap between 
the two beams gradually shrinks. Eventually the the ion 
two-stream instability grows to a non-linear stage, where 
large ion holes form and give rise to the ion density vari¬ 
ations seen in Eigure [2 

Simulations with shorter Debye length show 
instability-related density variations earlier in the 
wake, as can be seen in Eigure The large-scale 
structure of the wake is not sensitive to the change in 
Debye length, so the earlier appearance of the non-linear 
features of the two-stream instability is caused mainly 
by the faster growth associated with the higher plasma 
frequency at shorter Debye length; it takes a significant 
time for unstable perturbations to grow large enough to 
be visible when the Debye length is not infinitesimal. 
Extrapolating to very small Debye-lengths indicates 
that the ion beams become two-stream unstable when 
t ^ 3 ^m/c 3 , which corresponds to a distance of r^TbRu 
behind the moon for the particular solar wind speed 
used in the figures. 

Traditionally, Maxwellian beams are considered when 
studying the two-stream instability (e.g.. Ref. |23]), and 
cold ion beams are found to be unstable when their sep¬ 
aration is less than about twice the sound speed. How¬ 
ever, as illustrated in Eigure the ion beams in the sim¬ 
ulated wake are not Maxwellian, and have sharp edges 
associated with the initial void. To include enough par¬ 
ticles to resolve the fine features in velocity-space, large 
spatial bins of width dAoe are used when storing the dis¬ 
tributions for the following stability analysis, leading to 
a small amount of smoothing. 

The electrostatic ion stability at each location in 
the wake has been analyzed by calculating numerically 
the complex Eourier-transformed susceptibility of the 
(slightly further smoothed) ion distribution at that lo¬ 
cation together with Maxwellian electrons, in a ID lin¬ 
earized Vlasov-Poisson system, and setting it equal to 
zero. The ion-electron instability is suppressed in this 
calculation by artificially setting the imaginary part of 
the electron susceptibility to zero; this is appropriate 
because a Boltzmann electron treatment excludes elec¬ 
tron Landau damping or growth. Important deviations 
of the electron distribution from a Maxwellian will be 
discussed later, but are not considered here as atten¬ 
tion is restricted to the ion two-stream instability. The 
growth rate 7 of the most unstable mode at each location 
is shown in Eigure which clearly supports the interpre¬ 
tation of Eigureand the inferred instability onset. The 
position chosen for the example distribution in Eigure 
can be seen to lie close to the boundary between stable 
and unstable distributions. 

If the imaginary part of the electron susceptibility is in¬ 
cluded, weak ion-electron instability appears across most 
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FIG. 1. Ion density ni (in units of solar wind density rioo), normalized charge density ni —ne, and potential 0 of the wake 
simulated with Boltzmann electrons for Aoe = 0.02i?M (with the illustrative plasma flow speed Vsw = 25cs). 



FIG. 2. Ion distribution at three different times (and thus x) 
in the simulation shown in Figure^ 


of the wake. This instability is likely associated with the 
sharp edges of the ion beams, and appears to be the 
instability discussed by Farrell et al.l^ The Boltzmann- 
electron simulations of this section do not capture that 
instability, but even in the fully kinetic simulations in 
Section [IV| it does not appear to play a large role. 

Before moving to fully kinetic simulations, it is use¬ 
ful to examine the electron distributions that arise from 
the potentials of the Boltzmann-electron simulations. 
HutchinsorP^ showed that for a model potential evolution 
the electron distribution develops an unstable depression 
(a dimple) at a position-dependent velocity. Such a dim¬ 


ple is indeed seen in the electron distributions from the 
Boltzmann-electron simulations, examples of which are 
shown in Figure The dimple is more pronounced for 
more artificial mass ratios, leading Birch & ChapmaiP 
to treat the electron distribution as two electron beams 
in their stability analysis for = 20. For a realistic 

mass ratio the dimple only affects a very narrow velocity 
range of an otherwise largely Maxwellian electron distri¬ 
bution, so treating the distribution as two beams is not 
appropriate. Also seen in Figure is the initial distur¬ 
bance to the electron distribution from the (somewhat 
unphysical) sudden removal of electrons from the com¬ 
putational domain. The disturbance propagates out of 
the domain very quickly at realistic mass ratios, so any 
inaccuracy associated with the ID treatment of the pas¬ 
sage of the moon is unlikely to affect the parts of the 
wake of interest in the present work. 


IV. FULLY KINETIC SIMULATIONS 

Protons are 1836 times more massive than electrons, 
providing a lower bound for the true ion to electron mass 
ratio in a typical plasma. Simulations run at such 

mass ratios are computationally expensive, so a common 
approach is to use an unphysically small mass ratio in 
an attempt to capture the essential physics of the scale 
separation without paying the full computational cost. 
This was done in prior PIC simulations of the lunar 
wak^3IISl^ which used the highly unphysical mass ratio 
= 20. Later, Hut chinsorP^ showed that an unphys¬ 
ical mass ratio greatly enhances the degree to which the 
lunar wake electron distribution is unstable, and Hong 
et al.l^ showed that an unphysical mass ratio can signif¬ 
icantly affect ion-driven beam instabilities in ways that 
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FIG. 3. Ion density in the wake for Boltzmann-electron simulations at three different Aoe- 



FIG. 4. Ion distribution at a: = 82.5i?M and y = ^Rm in 
Figure with Gaussian counting errors for each bin shown 
in red. 


are not easily scaled to a physical mass ratio. Acknowl¬ 
edging these findings, the present work uses simulations 
at unphysical mass ratios only to illustrate phenomena 
observed in simulations using the physical mass ratio for 
hydrogen. 

At the true mass ratio and Aoe = 0.0025i?M5 treating 
both electrons and ions kinetically in ESPIC gives rise 
to the ion density shown in Figure Treating the elec¬ 
trons kinetically gives rise to instability much earlier than 
was seen in the simulations with Boltzmann electrons (cf. 
Figure]^. Non-linear disruption of the ion beams in this 
case occurs at t ~ which corresponds to a dis¬ 

tance of ^ 28i?M behind the moon (at the particular solar 
wind speed used for illustration). The nature of the den¬ 
sity variations is also different, with fewer, less periodic 
density enhancements gradually widening. 

For the same simulation, the evolution of the electron 
distribution is shown in Figure Very early in the wake, 
the only noticeable features are the initial void and the 
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FIG. 5. Maximum linear ion-instability growth rate in wake 
for simulation shown in Figure in units of the solar wind 
ion plasma frequency cjpioo- 

(faint) outward propagating initial disturbance from the 
sudden removal of electrons. Later in the wake a narrow 
dimple is seen, showing that the dimple persists in some 
form also when self-consistently considering electron in¬ 
stabilities (cf. Figure]^. Eventually, electron holes with 
large velocity extent form and gradually widen in spatial 
extent, corresponding to the density enhancements seen 
in Eigurej^ 

Because simulations at the true mass ratio have an 
extremely narrow dimple and are somewhat noisy, a sim- 
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FIG. 6. Electron distribution from potential evolution in 
Boltzmann-electron simulation, at time corresponding lo x — 
15i?M, for two different mass ratios. 



FIG. 7. Ion density ni and normalized charge density ni—rie 
of wake simulated with kinetic electrons for Aoe = 0.0025i?M. 
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FIG. 8. Electron distribution at three different times (and 
thus x) in the simulation shown in Figure]^ 
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FIG. 9. Departure of electron distribution from a Maxwellian 
for a simulation with = 459 and Aoe = 0.0025i?M. 


holes. The ion and electron distributions in the central 
region of the domain are shown in Figure for a time 
when the hole has grown to a large enough extent to sig¬ 
nificantly perturb the ion beams. After that time the 
ion beams are quickly disrupted (similarly to in prior 
simulation^^^^ , further widening the electron hole in ve¬ 
locity as the gap between the ion beams closes. A mech¬ 
anism has been proposed to explain the growth of such 


ulation using = 459 (but the same Debye length) is 
now used to illustrate the relevant electron phenomena. 
To better visualize the dimple and holes, a Maxwellian of 
the same density (and the solar wind electron tempera¬ 
ture) is subtracted from the electron distribution at each 
point in space. The result is shown in Figurefor a time 
before electron holes with large velocity-extent develop. 
Small hole-like perturbations are seen along the dimple 
in the electron distribution, and examining consecutive 
time-steps reveals that they convect outwards along the 
dimple in phase space (without significantly perturbing 
the ions). 

At a slightly later time, an almost stationary electron 
hole near ^ = 0 starts growing in velocity extent; the 
hole moves slowly {vy ^ — Cg) towards negative y, not ap¬ 
pearing to accelerate to follow the dimple like the other 



FIG. 10. Ion and electron distributions at a somewhat later 
time than in Figure 
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hole^, whereby the increasing local density as the hole 
moves downstream drives an expansion of its velocity- 
extent. 

Returning to the simulation at the true mass ratio, 
a careful examination of the evolution of the electron 
distribution reveals that the large electron holes seen in 
Figure arise in the same way as the one discussed for 
the slightly reduced mass ratio: small holes emerge from 
the dimple near ^ = 0, move slowly away from the dim¬ 
ple while growing in velocity-extent and remaining cen¬ 
tered at small Vy^ and then eventually expand to their 
full velocity-extent as they disrupt the ion beams; the 
holes then keep expanding in spatial extent as seen in 
Figure The first seed of a large hole emerges at a time 
corresponding to x ~ where the local electron den¬ 

sity is 0.21noo, and the holes grow until they disrupt the 
ions around x = 28Rm, at which point the local density 
is in the range 0.46noo to 0.56noo (depending on the hole 
position). For the same time-interval the total velocity- 
separation of the ion beams shrinks from ^5cs to ^Scg. 
That density increase and beam-separation decrease is 
sufficient for the proposed hole growth mechanisrrP^ to 
lead to ion-beam disruption, consistent with the simula¬ 
tions. 

Because the electron holes do not remain exactly at 
^ = 0 the ion beams are not quite balanced, as illustrated 
in the (extreme) example in Figure]^ Further, the finite 
velocities of the electron holes mean that they are likely 
closer to one of the beams in velocity than to the other, 
and those velocity-separations evolve as time passes. The 
situation in the simulations is thus not quite that of the 
symmetric cold beams considered in the proposed hole 
growth mechanisrrf^, but the theory nonetheless seems 
to capture the main features of the hole growth and beam 
disruption. Where in the wake the ion beams are dis¬ 
rupted thus depends on two main factors: the amount 
of electron hole growth from the increasing local density, 
and the decreasing hole size needed to disrupt the ion 
beams as their velocity-separation shrinks. 

Further simulations have been carried out to explore 
the sensitivity of the ion-beam disruption to changes in 
mass ratio and electron Debye length. In general, de¬ 
creasing the Debye length or the ion to electron mass 
ratio leads to earlier disruption of the ion beams. For 
example, changing the mass ratio to = 115 leads 

to the ion beams being disrupted at a time corresponding 
to X ~ 12 Rm5 when they are separated by ^S.Scg. It is 
presently not computationally feasible to do true mass 
ratio simulations with ESPIC at Debye lengths much 
smaller than the Aoe = 0.0025Rm of the simulations 
presented in this section, so the Aoe 10 ^Rm of the 
solar wind remains out of reach. However, the observed 
trend is for there to be more electron holes at shorter 
Debye lengths, appearing and disrupting the ion beams 
even earlier in the wake; at Aoe = 0.00125Rm disruption 
occurs at X ~ 20Rm- 


V. OBSERVATIONAL SIGNATURES 

The ARTEMIS missiorP has been observing the lu¬ 
nar wake for several years, and a numbe r of pote ntially 
relevant observations have been presentSolar 
wind variablitity and short Debye lengths in the wake 
make direct comparisons with the simulations difficult, 
but a qualitative discussion of possible observational sig¬ 
natures of electron holes and ion-beam disruption is given 
in this section. 

Small electron holes convecting along the dimple in 
phase space are one of the two main kinetic phenom¬ 
ena observed in the simulations; the holes arise because 
the dimple is unstabl^^, and may be present throughout 
most of the observed wake. Since the holes remain close 
to the dimple in phase space their propagation speeds 
are of order the electron thermal speed (except very close 
to the center of the wake), and the associated pertuba- 
tions extend ^ lOAoe ^ space (parallel to the magnetic 
field) and ^^e/ 20 e in potential. Such a hole may thus 
for instance extend 100 m along the magnetic field, be 
moving with a parallel speed of 10^ ™/s, and have a po¬ 
tential amplitude of 0.5 V. That would give an electric 
field signal of ^10™^/m with duration ^0.1 ms, which 
though weaker and faster than thos e of e lectron holes 
observed in the Earth’s plasma sheefP^J^ may still be 
detectable with ARTEMIS. High-frequency electrostatic 
fluctuations have been observed in the lunar wak^^, and 
could in principle be poorly resolved electron holes. 

The other main kinetic phenomenon in the simulations 
is slow-moving electron holes that grow large enough to 
disrupt the ion beams; their typical propagation speeds 
are of order the sound speed, their potential pertuba- 
tions are and there is no clear upper limit on 

their sizes. Such holes may be present beyond a cer¬ 
tain distance behind the moon, which depends on the 
(variable) solar wind parameters but is expected to fall 
within the range covered by ARTEMIS. Since the ion- 
disrupting electron holes are slower and perhaps more 
extended than the small fast holes, their temporal signa¬ 
tures should be of significantly longer durations. That 
could make them easier to detect and characterise, but 
may also make them more difficult to disentangle from 
other signal variations. One of the high-resolution ob¬ 
servations reported by Tao et al.^^ revealed broad-band 
electrostatic fluctuations at around the ion plasma fre¬ 
quency, for which electron holes were proposed as one 
possible cause (cf. Ref. [STJ. Though no well-defined soli¬ 
tary structures were identified, the fluctuations could be 
a sign of ion-disrupting holes, or perhaps of small holes 
spawned close enough to the center of the wake to be 
slow. 

Halekas et al.^^ presented wake-scale data for the same 
traversal studied by Tao et al. Ubiquitous ion-disrupting 
holes could in principle affect the ion distribution even at 
that scale, but the extent to which such holes could have 
developed only ^ 3.5 Rm downstream of the moon is not 
known from the simulations. The observed ion spectro- 





gram is dominated by the solar wind flow, but does show 
signs of counterstreaming ions. A higher-resolution seg¬ 
ment (Figure 4 of Ref. [3]) could be interpreted as showing 
the gap between the ion beams closing (the beams be¬ 
ing disrupted) in some regions, but solar wind variability 
makes direct interpretation difficultPl. That the regions 
with signs of ion-beam disruption lie close to the center 
of the wake is suggestive, however, since that is where the 
seeds of the ion-disrupting holes originate in the simula¬ 
tions. It is also where the high-resolution data showed 
possible signs of slow-moving electron holeJ^, and ap¬ 
pears to be the source of an outward-propagating electron 
beanJ^ electron holes generally increase the parallel elec¬ 
tron temperature, and deep ones can appear as beams. 

There is now an abundance of ARTEMIS data avail¬ 
able, as illustrated by recent worlP^ creating a 3D map of 
the typical large-scale structure of the lunar wake. Possi¬ 
ble signs of electron holes have been discussed in this sec¬ 
tion, but specific searches and detailed analysis of high- 
resolution data are needed to confirm the presence of such 
holes in the lunar wake. 


VI. CONCLUSIONS 

The present work highlights the importance of using a 
realistic ion to electron mass ratio in ID simulations of 
magnetized supersonic wakes behind non-magnetic ob¬ 
jects. Doing so has revealed a novel phenomenon where 
the velocity extents of electron holes grow from the small 
scale associated with electron-electron instability to one 
where the holes can significantly perturb the ions. As a 
result, the dimple in the electron distribution remains im¬ 
portant despite its narrow width at realistic mass ratios, 
acting as a seed of electron holes which eventually dis¬ 
rupt the counterstreaming ion beams in the wake much 
earlier than in the Boltzma nn-e lectron simulations. 

Prior kinetic simnlation d^ * ^^ * observed strong kinetic 
electron effects at a highly artificial mass ratio ("^i/me = 
20), but could not establish the importance of those 
effects at physical mass ratios because of the strong 
mass-ratio dependence of the dimpl^I^. The present re¬ 
sults demonstrate unequivocally the importance of ki¬ 
netic electron effects in electrostatic ID simulations, thus 
calling into question the validity of the hybrid simulations 
frequently used in attempts to explain experimental ob¬ 
servations. It is conceivable that other physical effects 
suppress the kinetic electron effects: for example, mag¬ 
netic perturbations have been shown to be more impor¬ 
tant for solar wind flow which is far from perpendicular 
to the magnetic fielcP^, and so could in principle inter¬ 
fere with the parallel electron and ion dynamics. Unless 
such a suppressing mechanism can be identified, however, 
accounting for kinetic electron effects on ion stability ap¬ 
pears to be essential. 

Because of the small size of the ion thermal gyroradius 
compared with the radius of the moon, the ID treat¬ 
ment is likely resonably accurate for the early large-scale 


structure of the wake. However, the Debye length is sig¬ 
nificantly smaller than even the electron thermal gyro- 
radius, so the accuracy of the ID non-linear instability 
evolution is less certain. In principle the electrostatic 
modes present in ID could couple to electromagnetic 
modes, possibly at oblique angles, giving the instabili¬ 
ties a multi-dimensional nature and introducing pitch- 
angle scattering, which could affect the electron dimple. 
Further, the growing electron holes may not be stable 
for the true ordering of length scales, possibly interfer¬ 
ing with the ability of electron holes to reach the size 
needed to disrupt the ion beams. Kinetic simulations of 
higher dimensionality are thus needed, though these will 
be computationally challenging given the need to use a 
realistic mass ratio and a short Debye length. 

Though the lunar wake has been used as the main ex¬ 
ample of possible applicability of the present simulations, 
they solve the general problem of supersonic flow past a 
non-magnetic object. As such they apply equally well to 
probes in magnetized plasmas and other large bodies in 
the solar wind, provided the main assumptions hold. The 
importance of electron kinetic effects at supersonic flow 
also raises the question of the extent to which similar phe¬ 
nomena can occur at slower magnetized plasma flow, but 
to address that question will require multi-dimensional 
simulations since the finite extent of the object in the 
flow direction and its boundary conditions then matter. 
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